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[ Runaway solutions can be avoided in fourth order gravity by a doubling 

I of the matter operator algebra with a symmetry constraint with respect to 

l> 

' the exchange of observable and hidden degrees of freedom together with the 

(N ■ 
O 



change in sign of the ghost and the dilaton fields. The theory is classically 

o : 

equivalent to Einstein gravity, while its non-unitary Newtonian limit is shown 
to lead to a sharp transition, around 10^^ proton masses, from the wavelike 
properties of microscopic particles to the classical behavior of macroscopic 



X 

H 

\ bodies, as well as to a trans-Planckian regularization of collapse singularities. 

A unified reading of ordinary and black hole entropy emerges as entangle- 
ment entropy with hidden degrees of freedom. The emergent picture gives a 

substantial agreement with B-H entropy and Hawking temperature. 
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I. INTRODUCTION 



A possible way out of the so called information loss paradox |T|,^ emerging from black 
hole physics consists in assuming a fundamental non-unitarity In fact it is natural to 

expect that the decoherence due to black hole formation and evaporation should give rise to 
a significant modification of the dynamical evolution laws: " For almost any initial quantum 
state, one would expect there to be a nonvanishing amplitude for black hole formation and 
evaporation to occur - at least at a highly microscopic (e.g., Planckian) scale - thereby 
giving rise to a nonvanishing probability for evolution from pure states to mixed states" . 
Although such an evolution is incompatible with a cherished principle of quantum theory, 
which postulates a unitary time evolution of a state vector in a Hilbert space, the crucial 
issue is to assess if it necessarily gives rise to a loss of quantum coherence or to violations 
of energy-momentum conservation so large as to be incompatible with ordinary laboratory 
physics [^|-^,||. Arguments for such violations were given, starting from the assumption that 
the effective evolution law governing laboratory physics has a Markovian character @J^. On 
the contrary one would expect that an effective evolution law modeling the process of black 
hole formation and evaporation, far from being local in time, should retain a long term 
"memory" In particular the basic idea of the non-Markovian models considered in 

Ref. is to have the given system interacting with a "hidden system" with "no energy of 
its own and therefore... not... available as either a net source or a sink of energy". 

On the other hand a mechanism for large entropy production in gravitational collapses 
should most naturally operate in the high curvature region, where one may expect new 
physics to emerge, while connecting it with the event horizon is somehow puzzling, as the 
physics on such a manifold has nothing peculiar for a free falling observer. Of course a 
quantitative model of Bekenstein-Hawking (B-H) entropy , along these lines, has to refer 
to the collapsed matter and, in order to do that, it has to include a mechanism for the 
elimination of the singularity. This does not mean that one can not identify the entropy 
carried by Hawking radiation as coming from the horizon within a local viewpoint: the 
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entropy growth outside the horizon, instead of being directly connected with an entropy 
produced by a strongly non-unitary dynamics in the region close to the classical singularity, 
is locally seen as a transformation of entanglement entropy into von Neumann entropy. In 
fact the relative character of the degrees of freedom involved in a given entanglement entropy 
is present even in flat space-times, where it can be exhibited explicitly |T^. Of course, in 



trying to pass from the region close to the horizon, where conventional quantum field theory 
in curved space-times is expected to work as a good approximation, to the region close to 
the classical singularity, we have to pay a price. In the absence of a full theory of quantum 
gravity, we have to rely on partially heuristic arguments and some guessing work, which we 
intend to show can be carried out by rather natural assumptions. 

In looking for a non-unitary theory avoiding the collapse singularity, we are going to start 
from higher derivative (HD) gravity. From a purely cosmological viewpoint it achieved great 
popularity since an inflationary solution was obtained without invoking phase transitions in 



the very early universe, from a field equation containing only geometric terms |TT|. More 
recently a renewed attention towards HD gravity was sparked by the appearance of HD 
gravitational terms in the low-energy effective action of string theory and in the holographic 
renormalization group, as well as by a growing interest in the study of brane worlds in HD 
gravity [|12|. However, although HD theories of gravity are natural generalizations of Einstein 
gravity, already on the classical level they are unstable for the presence of negative energy 
fields giving rise to runaway solutions. On the quantum level, as to unitarity, a possible 
optimistic conclusion is that "the S-matrix will be nearly unitary [0]. The crucial 

obstacle in trying to define HD gravity as a sound physical theory, namely the presence of 
ghosts, seems in fact to be a strong indication, on one side, of non-unitarity and, on the 
other, of a possible mechanism for avoiding singularities, thanks to short range repulsive 
terms. 

Here a specific non-unitary realization of HD gravity is shown to be compatible with the 
wavelike properties of microscopic particles, as well as with the assumption of a gravity- 



induced emergence of classicality |15H23[| , and seems to give strong indications for the elim- 



ination of singularities on a trans-Planckian scale. Parenthetically we are encouraged in our 



extrapolations by the success of inflationary models, implicitly referring to these scales pi 
The present setting suggests that B-H entropy may be identified with the von Neumann en- 
tropy of the collapsed matter, or equivalently with the entanglement entropy between matter 
and hidden degrees of freedom, both close to the smoothed singularity. This viewpoint is 
corroborated by the attractive features of the Newtonian limit of the model. In this attempt 
of evaluating relevant physical quantities by an incomplete theory, and in particular by its 
Newtonian limit, we are encouraged by well-known precedents, like the amazing quantitative 
agreement between the analysis of John Mitchell in 1784 and the modern notion of a black 
hole. 

A bonus of the present non-unitary model is the possibility of a unified notion, as von 
Neumann entropy, both for B-H and ordinary thermodynamic entropy of closed systems. 
This is not irrelevant, as, " ...in order to gain a better understanding of the degrees of freedom 
responsible for black hole entropy, it will be necessary to achieve a deeper understanding of 
the notion of entropy itself. Even in fiat space-time, there is far from universal agreement 
as to the meaning of entropy - particularly in quantum theory - and as to the nature of the 
second law of thermodynamics" [^. 

II. STABLE FOURTH ORDER GRAVITY 



Long ago deWitt [|l^ and Stelle p5[ analyzed the improved ultraviolet behavior of HD 



gravity theories stemming from cancellations that are "analogous to the Pauli-Villars regu- 



larization of other field theories" . These cancellations are precisely due to the presence 
of negative energy fields, which in their turn are the source of instability: energy can flow 
from negative energy degrees of freedom to positive energy ones and one can have runaway 
solutions. 



In Ref. ||T4[ a remedy for the ghost problem, leading to a non-unitary theory, was sug- 
gested by a suitable redefinition of the euclidean path integral. In this paper we mean to 
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propose an approach directly in real space-time, thus avoiding analytic continuation, which 
amounts to a very tricky operation outside the realm of a fixed fiat geometry. Like in Ref. 



14 1 , classical instability is cured at the expense of unitarity. Before treating the physically 



relevant case, we consider first a simpler fourth order theory for a scalar field 0, which has 



the same ghostly behavior as HD gravity WM. Its action is 



S 



d X 



(1) 



with the inclusion of a matter action S^at and an interaction with matter, where ip^ip is a 
shorthand notation for a quadratic scalar expression in matter operators. Defining 



□0 



(2) 



the action can be rewritten as 



S [01, 02, Ip] = Smat 1p] 
+ I d\ 



^0in0i - (□ - /x^) 02 - A 



a 



+ -^V (02 



(3) 



The action of 02 has the wrong sign, which classically means that the energy of the 02 field 
is negative. If there were no interaction terms, this negative energy wouldn't matter because 
each of the fields would live in its own world and the positive and negative energy worlds 
would not communicate with each other. However, if there are interaction terms, like 0^ 
or '0^'005 energy can fiow from negative to positive energy degrees of freedom, and one can 
have runaway solutions, with the positive energy of 0i and the negative energy of 02 both 



increasing exponentially in time 14 



This toy model shares with HD gravity theories some of the mentioned improvements on 
the ultraviolet behavior. In fact there is a complete cancellation of all infinities coming from 
the '0^'00 interaction and corresponding to self-energy and vertex graphs I^G], owing to the 
difference in sign between 0i and 02 propagators. A key feature of the non-interacting theory 
(A = a = 0), making it classically viable, can be considered to be its symmetry under the 
transformation 02 — > —02) by which symmetrical initial conditions, i.e. with 02 = 02 = 0, 



produce symmetrical solutions, thus in particular avoiding the runaway ones. We are going 
now to extend this symmetry to the interacting theory. If one symmetrizes the action as 
it is, this eliminates the direct interaction between the ghost field and the matter altogether 
and then the mentioned cancellations. A possible procedure to get a symmetric action 



while keeping Pauli-Villars-like cancellations is suggested by previous attempts and by 
the information loss paradox both pointing to a non-unitary theory, where tracing out 
hidden degrees of freedom results in general in mixed states. In particular the most natural 
way to make the hidden degrees of freedom "not... available as either a net source or a sink 
of energy" is to constraint them to be an exact copy of the observable ones. Accordingly 
we introduce a (meta-)matter algebra that is the product of two equivalent copies of the 
observable matter algebra, respectively generated by the and ip^ip operators, and a 

symmetrized action 

Ssym=^\^S [4)i,(f)2,'ip\lp] + S 01, -02,^^V^ } , (4) 

which is invariant under the symmetry transformation 

h — ^ -02, -ip — ^ V", — ^ • (5) 
This duplication is formally analogous to what is done in thermo-field dynamics p7 



where in particular it can be used to describe the irreversible evolution of open systems p8 



If the symmetry constraint is imposed on quantum states, i.e. the state space is restricted 
to those states |^) that are generated from the vacuum by symmetrical operators, then 

This implies that, as usual with constrained theories, allowed states do not give a faithful 
representation of the original algebra, which is then larger than the observable algebra. In 
particular the constrained state space cannot distinguish between F [ip\ip\ and F 
by which the operators are referred to hidden degrees of freedom, according to a standard 
terminology for non-unitary models |^ , while only the ip operators represent matter degrees 
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of freedom. On a classical level the constraint implies that and ip are to be identified 
while the 02 field vanishes and, as a consequence, the classical constrained action is that of 
an ordinary second order scalar theory interacting with matter: 



ci 



d X 



)in0i - A 



/i 



a 



/i 



(7) 



Consider then the classical action for a fourth order theory of gravity including matter 



m 



S = Sg [g^^u] + Smat [Qi^u, tp\ ^] 



-9 



1 



:7^ 



d'^x\^~gL 



mat J 



(8) 



16nG 

where Lmat denotes the matter Lagrangian density in a generally invariant form. In terms 
of the contravariant metric density \/'i2'KGh^^ = \/—gg^^ — rj^^ the Newtonian limit of the 
static field gives 

1 1 e"^"'' 4 e-'^^r 



00 



o o , (9) 

I 6 r 6 r 

where /xq = [327rG(3/3 — /i2 = [IGttG'q;]^"'^/^ while a complete analysis for the 

whole metric can be found in Ref. From Stelle's linearized analysis, the first term in 
Eq. (|^) corresponds to the usual massless graviton, the second one to a massive scalar and 
the third one to a negative energy spin-two field. In fact, in analogy with Eq. one can 
introduce an explicit transformation from the initial field g^y appearing in the fourth order 
form of the action to a new metric tensor g^^, a massless scalar field x dilatonically coupled 
to the metric and a spin-two massive field (f)^y, this transformation leading to the second 
order form [Q. To be specific, following Ref. (see Eq. (6.9) apart from the matter 
term), the action (P) can be rewritten as the sum of the Einstein-Hilbert action Seh for g^^p, 
an action Sgh for the traceless symmetric ghost field 0^,^ and the scalar field x coupled to 
the metric g^^, and a matter action Smat^ where g^y is expressed in terms of g^y, (p^y and x 
(replacing g^y by e^g^y in Eq. (4.12) in Ref. Q): 

S [g^y,(j)^y,X,'^\i^] 

= Seh [gfiu] + Sgh [(jfiu, <Pfj.u, x] + Smat [gfj.u{g x),^^V']• (10) 
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In Sgh above the quadratic part in 0^,^ has the wrong sign just as for 02 in Eq. (|^). As 
before, a simple way to get rid of classical instability would be to symmetrize the action with 
respect to the transformation (p^^, —(pfiu and to introduce the symmetry constraint with 
respect to this transformation. This however would eliminate the corresponding repulsive 
term in Eq. (^, which is a possible candidate in avoiding the singularity in gravitational 
collapse. Once one accepts non-unitarity, it is rather natural to assume that one can cure 
the instability, while keeping the short-range repulsive term, by introducing hidden degrees 
of freedom as above, i.e. from a quantum viewpoint to accept that the operator algebra 
involved in defining the dynamics is larger than the observable algebra. To be specific, once 
again, we double the matter algebra by taking a meta-matter algebra which is the product 
of two copies of the observable matter algebra, respectively generated by the operators , ip 
and We then define the symmetrized action [plU 



(11) 



which is symmetric under the transformation 



-X, 4^ — ^ ^, — ' i^, 9t^u 



(12) 



Like above, if the state space is restricted to those states that are generated from the 
vacuum by symmetrical operators, then 



(13) 



Like for the previous toy model, the constrained state space does not distinguish between 
F '0] and F '0^, -0 , by which the ip operators are referred to hidden degrees of freedom, 
while only the ip operators represent observable matter. On a classical level the constraint 
implies that ip and ip are to be identified, while the (p^i, and x fields vanish and, as a 
consequence, the classical constrained action is that of ordinary matter coupled to ordinary 
gravity: 



(14) 
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as Sgh[g^..,^A = (Eq. (6.9) in Ref. 0) and 0, 0) = g^, (Eq. (4.12) in Ref. 

pO| with e^g^y replacing g^^). While this modification of fourth order gravity is expected 
to affect its ultraviolet behavior, still it does not worsen it for one-loop meta-matter to 
meta-matter amplitudes, at variance with the trivial symmetrization of the original theory. 
It should also be remarked that one could limit symmetrization to the ghost field only, 
without involving the scalar field, especially if one were concerned with the cosmological 
implications of keeping the dilatonic scalar field in the classical action. 

A final remark is in order as to the possible rereading of the present model as a bimetric 
HD theory where two worlds interact only by means of a coupling between the metrics 
(of the fourth order formalism) ||32|. In fact the model can be defined by replacing g^^, 
—<f)fj,u, and —X in the second term of the right hand side of Eq. by three independent 
fields and adding an interaction, including Lagrange multipliers, leading to the necessary 
identifications. However we do not commit ourselves to the prevailing view pointing to 
underlying higher dimensional theories, even though a natural setting could appear to be an 



extension of the Randall-Sundrum model with two positive tension fiat branes separated 
by one intermediate negative tension fiat brane , where ip and ip meta-matters reside on 
distinct positive tension branes. 



III. NEWTONIAN LIMIT AND GRAVITATIONAL LOCALIZATION 

Of course the elimination of classical runaway solutions is only a first step in assessing the 
consistency of the ensuing non-unitary theory. A further natural step consists in studying 
its main implications for ordinary laboratory physics. In order to do that, consider the 
Newtonian limit of such a theory with non-relativistic meta-matter and instantaneous action 
at a distance interactions. Looking at the signs in Eq. (|11|), we see the following. The 
interactions due to the massless graviton field cjf^i, are always attractive, whereas those due to 
the scalar field x are attractive but for the ones between observable and hidden meta-matter; 
finally those due to the massive field are repulsive within observable and within hidden 



9 



meta- matter, due to its ghostly character (see the sign of the third term in Eq. (^)), and 
are otherwise attractive, since the ghostly character is offset by the difference in sign in its 
coupling with observable and hidden meta-matter. The corresponding (meta-)Hamiltonian 
operator is then 

2^ J \x-y\ \ 3 3 

G f '^]{x)ipi{x)'ipl{y)ipkiy) ( 1 I I 4 II 

4^ y \x-y\ V3 3 

4^i \x-y\ V3 3 J 

acting on the product ® of the Fock spaces of the (non-relativistic counterparts of 
the) ip and ip operators. Here two couples of non-relativistic meta-matter operators 'ippipj 
and appear for every particle species and spin component, while nij is the mass 

of the j-th particle species and Hq is the matter Hamiltonian in the absence of gravity. 
The ipj operator obeys the same statistics as the corresponding operators ipj, while 
_ = = 0. Though never appearing in our formulae, the electromagnetic potential 

in the Coulomb gauge should be included in the original degrees of freedom, even though, 
in the non-relativistic setting, it is not involved in the gravitational interaction. 



With reference to Eq. (0), observe that the action at a distance counterpart of the 
field-theoretic cancellations mentioned above is the possibility of avoiding normal ordering 
in the last two terms. It would correspond, in fact, to the subtraction of the finite oper- 
ator G(/io — 4yU2) X]j "^j / dxipj{x)iljj{x)/ 12 and its hidden correspondent, which in a fixed 
particle-number space correspond to irrelevant finite constants. 

To be specific, the meta-particle state space S is the subspace of F^ (g) F^ including the 
meta-states obtained from the vacuum 1 10)) = |0)^® |0)^ by applying operators built in terms 
of the products ipj{x)ijjj{y) and symmetrical with respect to the interchange ip"^ ^ ip"^, which, 
then, have the same number of ip and ip meta-particles of each species. As the observable 
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algebra is identified with the if) operator algebra, expectation values can be evaluated by 
preliminarily tracing out the if) operators. In particular, for instance, the most general 
meta-state corresponding to one particle states is represented by 



11/))= / dx I rfy/(x,y)^](x)V^J(y)||0)), f{x,y) = f{y,x). 



(16) 



This is a consistent definition since Hq generates a group of (unitary) endomorphisms of S. 
A pure n-particle state, represented in the traditional setting by 



g) = j rf"x^(xi,X2,...,a:„)V^]^(xi)V^j2(x2)...V^],^(x„) |0) 



(17) 



is represented in S by the only meta-state that, by tracing out ip operators, gives the state 
\g) {g\, with \g) as in Eq. (|T3), namely by 

11^®^)) oc / d''xd'^yg{xi,...,Xn)g{yi,...,yn)i']^{xi)...%l)]^{xn)^^{yi)...^^^^ (18) 



It should be remarked that, when our initial knowledge of the system state is char- 
acterized by a density matrix, there is no unique prescription to associate it with a pure 
meta-state. In such a case one has to consider the possibility of using mixed meta-states to 
encode our incomplete knowledge. 

Considering, for notational simplicity, particles of one and the same species, the time 
derivative of the matter canonical momentum in a space region VL in the Heisenberg picture 
reads 



G=0 



+ FG = --[-fn,Ho[^lj\ 



+ -m^ I dx^\x)i,{x)V. I rfy^^MiM ( 1 _ le-^ol-j/l + Ig-M^k-.l 
2 Jn Jr3 F-y| V 3 3 

2 Jn Jr3 \x-y\ \ S 3 



(19) 



The expectation of the gravitational force can be written as 



2 \Jq Jn F — 1/1 V 3 3 

^ [ dx^HxWxW. ! d.^i^ ( 1 _ le-ol^'-l + i 



2 \ic iR3\f7 V 3 3 

+ / / dx^\x)i.{x)V. ! rf^^^M^ f 1 + le-.oN-.l _ le"-!-^!') \ 
2 \in Jo |x -1/1 V 3 3 7/ 

+ ^^2 / f rf^^t(^)^(3,)v, / rf^^lMM ( 1 + ie-Mo|.-.l _ i^-,.y^^y\\ \ ^ (20) 
2 \Jn JR3\n \x-y\ \ S 3 // 

where, on allowed states, the first term vanishes for the antisymmetry of the ker- 
nel Vx [(l — e~^'''^~^'/3 + 4e~'^^'^~^'/3) / |x — ?/|] and the symmetry constraint on the 
state, while the third one vanishes, as it should be for self-gravitating matter, just 
as a consequence of the antisymmetry of the corresponding kernel. As is usual 
with the evaluation of forces between macroscopic bodies, we can then approximate 
(^"^ {x)il){x)'ip\y)'ip{y)^ and (^il)\x)il){x)il)\y)il){y)') respectively by l^ip'^ {x)il){x)') {^\y)'ip{y)^ 
and i^ip'^ {x)il){x)') {ip\y)ip{y)) , as x G and y G R^\Vt. Finally, as {i'\y)'ip{y)^ = 
{ip^{y)ip{y)), we get 

(Fo) ^Gm" ! dx{ilj\x)^{x))W^ ! c/yi^^MlM, (21) 
^ ' Jn JR3\n F ~ y\ 

namely the classical aspects of the interaction are the same as for the traditional Newton 
interaction, consistently with the classical equivalence of the original theory to Einstein 



gravity pT| 



Although we are using the general Newtonian limit (0), it is worthwhile to remark that 
we are mainly interested to two opposite specialized limits. 

The ordinary Newtonian limit, for ordinary laboratory physics, corresponds to taking 

1^0,1^2 oo, if fiQ^ and /i^^ are assumed, as usual, of the order of the Planck length, in 

which case the meta-Hamiltonian Hq can be rewritten in the form 

~, ~ G sr^ f i)\{x)i)Ax)^Uy)i)k{y) 

^ ■ , J \x y\ 

where H[%p\%p\ and H[^\ip\ respectively include the halved (normal ordered) Newton in- 
teraction within observable and hidden meta-matter. In this form we have a well defined 
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non-unitary model of Newtonian gravity without any free parameter. Tracing out the ip op- 
erators from the meta-state evolving according to the unitary meta-dynamics generated by 
Hg results in a non-Markov non-unitary physical dynamics for the ordinary matter algebra 
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The trans-Planckian Newtonian limit concerns the use we are going to make of the model 
with reference to gravitational collapse, where the model replaces the classical singularity 
with a trans-Planckian structure. To this end we consider the opposite limit fio,fi2 —>■ 0, 



leading to a Hamiltonian Hq as in (p2|) with H and G/2 respectively replaced by Hq and 
G. The rationale for the use of this bold extension of the Newtonian limit outside its 
typical applicability range, though within a merely heuristic approach, resides in part in the 
soundness of its physical consequences, as shown in the following. 

A general new feature of the model with respect to the usual inclusion of Newtonian grav- 
ity in QM is the localization due to the presence of an effective self-interaction. Consider in 
fact in the traditional setting a physical body in a given quantum state whose wave function 

the product of the wave function of the center of mass and of an 
internal wave function. In particular "^cm can be chosen, for simplicity, in such a way that 
the corresponding meta-wave function ^tot = ^^CMi^)"^ iNT{xi — Xj)'ifcM{y)'^ iNTiUi — yj) 
can be rewritten as: 

^TOT = ^Cm( )^INTiX - Y)'^jNTiXi - Xj)^iNTiyi " %), (23) 

where yi,Y denote the hidden correspondents of Xi,X. As to ^uyxiX — Y), we choose it 
as the ground state of the relative motion of the two interpenetrating meta-bodies, which 
is formally equivalent to the plasma oscillations of two opposite charge distributions. The 
corresponding potential energy, if the body is spherically symmetric and not too far from 
being a homogeneous distribution of radius H and mass M, has the form ^GM'^f {\X — Y\), 
where 

— 1/r for r > 2H 
fir)={ ^, " ^, (24) 

iar^/S'^ for r ^ H 
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with ^ = 1/2, 1 respectively for the ordinary and the trans-Planckian hmit, and a ~ 10° a 
dimensionless constant. We are interested here to the case of small relative displacements. 
The relative ground state is represented by 

^iMX -Y) = (A^tt)"'/" e^^; A = {2n''Eya^GM^y/\ (25) 

Then, if we choose "^cAii^) oc exp [— X^/A^], we get 

^o(X, Y) = ^cm{X)^cm{Y) = ^int{X - Y)^jnt{X + Y). (26) 

In particular for body densities ~ lO'^^nip/cm^ , where rrip denotes the proton mass, A ~ 
{nip/ My ^'^ cm, which shows that the small displacement approximation is acceptable already 
for M ~ lO^^TUp, when A ~ 10^^ cm, whereas the body dimensions are ~ 10 ^cm p5| . 



Another simple case corresponds to masses lower than lO^'^mp, where the two meta- 
bodies can be approximated as point particles and their ground state wave function, in the 
ordinary Newtonian limit, is 

*(X - r) oc e-l^-^l/^ a = Ah^C^G-Hr^ ~ 10'^ {M/mp)-^ cm, (27) 

by which gravitational localization, consistently with recent experiments, can be ignored for 



all practical purposes even for particles much larger than fullerene ||36|j37|| . The ensuing sit- 
uation corresponds then to a rather sharp localization mass threshold Mt ~ h^^^G'^^^^ p^/^'^ , 
which is very robust with respect to mass density variation. 

It is easily seen that the present framework actually is compatible with the way terrestrial 
gravity appears in QM. A crucial experiment, dating back to 1975, exhibits in fact in a 
striking manner how terrestrial gravity enters the Schrodinger equation in the usual way, i.e. 
just as a Coulomb external field . To this end the calculation of the average gravitational 
force acting over a lump performed above does not suffice since it can explain only e.g. the 
free fall of a microscopic particle by means of classical equations (Ehrenfest theorem) where 
h does not appear. 
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Consider the problem of a large, for simplicity spherically symmetric, massive body (the 
Earth) in some irrelevant internal state in interaction with an external microscopic particle. 
Define the meta-Hamiltonian of the Earth-particle system as 

H=^Y.^%-GM^}m-R.\,-^Y.^l-Gn.M ^ ^ (28) 

where M and m are respectively the mass of the Earth and of the particle, it!i, R2 and xi, X2 
respectively the center of meta-mass coordinates of the two Earth and particle copies. 

Let's start with a meta-state of the meta-Earth system corresponding to the fundamental 
(or a not too highly excited ) one with respect to the relative motion of the two copies and 
choosing the initial CM meta-state of the bound system of the two copies just as above. The 
localization length is in this case of the order A^"''*'* ~ 10~^^ cm. Having in mind that the 
particle is described by a wave packet whose size a is in any case much larger than A^"''*'^, we 
can approximate the squared modulus of Earth's meta-wave-function by a product of delta 
functions 5^ (r) 5^ (R), where r, R respectively denote the internal and the CM coordinates 
of the meta-Earth bound system. As a consequence Xi — Rj in the Newton potential can 
be replaced by Xi. Of course, since the spreading time of the Earth's CM wave function 
over a region of the size a > 10"^*^ cm is given by aA^"''^^^M/h > 10^^ s, the approximation 
is justified in any physically relevant situation, and actually even much better than what 
appears from this analysis, as we are ignoring the spreading of the particle wave function. As 
a result the gravitational interaction enters in the particle dynamics simply by the presence 
of the usual external Newton potential. 

IV. EVOLUTION FROM PURE TO MIXED STATES 

It should be stressed that, while in the ensuing dynamics the constraint on the hidden 
degrees of freedom to have the same average energy as the observable matter avoids them 
to be "available as either a net source or a sink of energy", only the meta-Hamiltonian is 
strictly conserved. If we include in the physical energy the usual Newtonian interaction 
between observable degrees of freedom, the physical energy operator 
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Hph[ip\ij] = Ho[ip\tp] - — > mjUik I dxdy — . . 

2 ^-^ J F — y| 



(29) 



is not the generator of time evolution. To be specific, in the ordinary Newtonian hmit, the 
generator of the meta-dynamics can be written 



(30) 



^'^^mjUik I dxdy 



2iljUx)ijj{x)'iljl{y)iJk{y) 



\x-y\ 



+ — ^^m^mfc / dxdy 



■ '^j{x)i^j{x)'ipl{y)'^k{y) ■■ + : '4'j{x)ijj{x)'ipl{y)^k{y) ■■ 

\x - y\ 

from which we see that He and Hph[ip'^ + Hphlip^ip] in general are different only due 
to correlations. The two sums above have approximately equal expectations and fluctuate 
around the classical gravitational energy. On one side these energy fluctuations have to be 
present in any model leading to dynamical wave function localization, which in itself requires 
a certain injection of energy |^9|]. On the other hand these fluctuations, though irrelevant 
on a macroscopic scale, are precisely what can lead to thermodynamical equilibrium in a 
closed system if thermodynamic entropy is identified with von Neumann entropy [0. In fact, 
due to the interaction with the hidden degrees of freedom, a pure eigenstate of the ordinary 
energy Hph is expected to evolve into a microcanonical ensemble. 

As a simple example showing how a pure state can evolve into a mixed one, consider 
a free spherically symmetric body of ordinary matter above localization threshold, initially 
described by a gaussian wave packet, whose size is chosen as above in such a way that the 
particle-copy system is in its ground state, thus recovering the meta- wave-function (|23|) \\iO\\ . 

The factor depending on the center of meta-mass of the ip and ip meta-bodies in \E'o(-^, Y) 
(p6|), for M > lO^^rUp, spreads in time as usual for a body of mass 2M, so that after a time 
t, the meta-wavefunction becomes 



^t(X,y) ocexp 



\X -Y\ 
2A2 



exp 



- \X + Yf /A 



^ ^-ao\X-Y\' ^-at\X+Y\ 



(31) 



In order that this be compatible with the assumption that gravity continuously forces lo- 



calization 1 15 - 23 1 , the spreading of the physical state must be the outcome of the entropy 
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growth. This initially vanishes, as the initial meta-wavefunction ( ]26|) is unentangled and 
then the physical state, obtained by tracing out Y, is pure. If one evaluates the physi- 
cal state pt{X,X') = f dY'^t{X,Y)'^l{X',Y), one finds that the space probability density 
reads 



PtiX,X) 



8ao{at + at) 
n^at + at + 2ao) 



1 3/2 



exp 



Sapjat + at) 
(at + at + 2ao 



-X' 



cx: exp 



A4 + 2^2^VM2' 



(32) 



The spreading is extremely slow, as its typical time, for bodies of density ~ 10 nip/ cm 



is ~ lO^sec independently from the mass, as can be checked by means of Eqs. ( |52| ) and 
(^). If it is due to entropy growth only, rather than to the spreading of the wave function, 
the entropy St is expected to depend approximately on the ratio between the final and the 
initial space volumes roughly occupied by the two Gaussian densities, according to 



St ~ Kb- In 



at + at + 2ao 



(33) 



2{at + at) 

at least for large enough times. (Linear momentum probability density does not depend on 
time.) Of course this corresponds to the approximation of the mixed state by means of an 
ensemble of equiprobable localized states, which is legitimate if turns out to be large 
enough. In order to evaluate the entropy of the state represented by pt{X, X') and to check 
Eq. (|33D, we use the possibility, in this approximation, of linking the entropy 



St = -KBTr [ptlnpt] = KBlnN 



(34) 



with the purity 



Tr [p?] = i/AT; p2(X, X') = I dX"ptiX, X")pt{X\ X'). 



(35) 



By an explicit computation we get 
Tr[pl] = [ dXpl{X,X) = 



[4:ao{at + at)Y 



{2atat + QataQ + Qata^ + ^a^)^ — 4 [at — a^f {at — a^f 



and, for large times, one can keep just the leading term in aj/ao, that is 



3/2 

(36) 
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which, by using Eqs. (ISIJS^), gives 



S, . -k' .„ f ^) ^ A-4 „> f ^l±i^!^V^l! I (38) 



2 V 2ao y 2 V 

which differs from the leading term in Eq.(^3l) by an irrelevant quantity (3/2)i^'B In 2. This 
validates our view of the free motion of a macroscopic body, at variance with the rather 
unphysical stationary localized states of the Schrodinger-Newton (S-N) model, whose initial 
linear momentum uncertainty does not give rise to a spreading of the probability density 



41|^5[. More generally, while that non-linear generalization of QM was considered to be 
a reasonable mean field approximation of an unspecified theory, by its unitarity it can not 
model any fundamental gravitational decoherence. It is remarkable that the S-N model 
can be actually obtained as the N ^ oo limit of the color generalization of the present 



Newtonian limit [46|. 



It should be stressed that the notion of coarse graining entropy, often taken as the starting 
point in dealing with the quantum foundations of the second law of thermodynamics ||^, can 
be easily connected with the present approach. Consider, for simplicity, a non-degenerate 
physical state 

pph = ^Pj \j) 01 , Pj e R, Pj =Pk^j = k. (39) 
j 

The most general pure meta-state vector giving rise to pph is 

j 

where \ j) \ j) denotes the tensor product of two corresponding vectors in the two Fock spaces 
and the (pj G [0,27r[ are arbitrary real parameters. The indistinguishability of the corre- 
sponding meta-states, due to the restriction of the physical algebra, induces in the meta- 
state space an unambiguous coarse graining, at variance with the rather vague one in the 
traditional approaches. To be specific, it is natural to introduce the macro-meta-state 



p^^^ I \{'^m).m\ = Y.p^\j)\j) (41) 

corresponding to the equiprobability of the micro-meta-states ((^11 • The correspond- 

ing coarse graining entropy is 

ScG = -KbTt [pcG In Pcg] = -Kb ^ Pj Inpj , (42) 

i 

which coincides with the von Neumann entropy of the physical state pph- 

Vice versa, if we assume that a specific pure meta-state \\^)) is given, the Schmidt 
decomposition theorem allows us to write it in terms of orthonormal vectors as 

ll*))=Ev^l^')l^")' (43) 

with the Pj positive, for simplicity distinct, real numbers. By the symmetry constraint on 
the meta-state space one can choose the relative phases in such a way that \j) and 
can be taken as corresponding vectors in the two Fock spaces, thus reproducing H^E'))^ in 
eq. (JD) for y9 = 0. Although this amounts to the knowledge of a definite microstate, the 
entropy of the corresponding physical state pph is non- vanishing and coincides with the 
coarse graining entropy of the corresponding macrostate pcc- This shows the objective and 
non-conventional character of the notion of entropy in the present approach, since it does 
not depend on a subjective characterization based on the notion of a macroscopic observer 



V. WAVE FUNCTION REDUCTION 



In an interaction representation of the ordinary Newtonian limit, where the free meta- 
Hamiltonian is H[^\ip\ + Hl^p^if)], the time evolution of an initially untangled meta-state 
$(0)')') is represented by 



Texp 
U{t) $(0) 



— Gm^ I dt I dxdy 



\x - y\ 

uit)\m)^®\m)^- 



m 



(44) 
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Then, by a Stratonovich- Hubbard transformation [Q, we can rewrite U{t) as 

U{t) = j D [v?!, (p2] exp ^ [ dtdx [v^iV^i - 
Texp 
Texp 



-i^-V2TTG I dtdx[(fi{x,t) — (p2{x,t)]ilj'^{x,t)ijj{x,t) 



namely as a functional integral over two auxiliary real scalar fields ipi and ip2- 
The physical state corresponding to the meta-state (Q) is given by 



k 

and, by using Eq. (^Sj), we can write 

^(t))) = j ^ bi' <^2] exp ^ y dtdx [<^iVVi - <^2VV2] 



(45) 



(46) 



^ Texp 
Texp 



—i—-V2TiG / dtdx [v9i(x, t) — (^2(2;, t)] '?/'^(x, t)^/'(x, t) 
h 



i^-y/27cG I dtdx[(fi{x,t) + (f2{x,t)]ilj\x,t)ilj{x,t) 



mk- (47) 



Then the final expression for the physical state at time t is given by 

pphit) = I D[ipi,ip2, V^'i, '^'2] exp ^ / dtdx [v^i VVi - '^2^'^ ^2 - V^iVVi + V52VV2] 



($(0)1 r-^exp 



%—-\J1'kG \ dtdx \ip\ — (/Pol 
n 



Texp 
Texp 



— z- 



.mc 



■\/2ttG j dtdx [ifi — <~p2] i^^ip 



-i^-V2TTG I dtdx [<y9i + 992] 'i/'V 



1^(0))^ 
l*(0))^ 



($(0)1 T-^exp 



i^V2^ I dtdx[ifi[ + ip2]ip^ip 



(48) 



where, due to the constraint on the meta-state space, ip operators were replaced by 
operators, and the meta-state vector |$(0))^ by |$(0))^. This expression can even be taken 
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as an independent equivalent definition of tlie non-unitary dynamics, free from any reference 
to the extended algebra including unobservable degrees of freedom. 

Consider an initial linear, for simplicity orthogonal, superposition of localized states of 
a macroscopic body, existing, as shown above, as pure states corresponding to unentangled 
bound meta-states for bodies of ordinary density and a mass M higher than ~ lO^^mp 

N 



1 ^ 



(49) 



where \z) represents a localized state centered in z. We consider the localized states as 
approximate eigenstates of the particle density operator, i.e. ip^x, t)ip{x, t) \z) ~ n{x—z) \z), 
where time dependence is irrelevant, consistently with these states being stationary both in 
the gravity-free and in the interacting Schrodinger pictures apart from a slow spreading, 
which, as shown below, is much slower than the computed time for wave function reduction. 
According to Eq. (^), the density matrix elements are then given by 

{Zh\ PPh{t) \Zk) 

D [ipi, v?2, V^i, ^'2] exp ^ j dtdx [v^iVVi - V52VV2 - V^'iVVi + V52VV2] 

^ ^ V r 
-j^ exp —i^-V2nG J dtdx [[(fi — (P2] n{x — Zj) — [(p[ — (p'2] n{x — z. 

y/27rG J dtdx [[(pi + ^92] n{x - Zh) - [lp[ + (^'2] n{x - Zk)] 



exp 



.mc 



n 



(50) 



and, after integrating out the scalar fields. 



N 



{zh\pPh{t) \zk) = ]^X]exp-Gm^t I dxdy 



n{x - Zj)n{y - Zh) n{x - Zj)n{y - Zk) 



F - y\ 



\x-y\ 



(51) 



which shows that, while diagonal elements are given by {zh\pph{t) \zh) = 1/N, the coher- 
ences, under reasonable assumptions on the linear superposition in Eq. p9|) of a large 
number of localized states, approximately vanish, due to the random phases in Eq. (^TJ). 
This makes the state pph{t), for times t >Tg ^ 10^'' (M/mp)~^/^ sec, which are consistently 
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short with respect to the time of the entropic spreading ~ 10^ sec, equivalent to an ensemble 
of localized states: 

1 ^ 
i=i 

It is worthwhile to remark that the extremely short localization time of a macroscopic body 
may make its unlocalized states unobservable for all practical purposes. The above analysis 
is also supported by numerical evidence independently from the particular assumptions 
made here on the initial unlocalized state (^91) |50| . In such a way one gets a gravity- 



induced dynamical reduction of the wave function, which up to now was assumed to follow, 
possibly, from a future theory of quantum gravity [0. It is worthwhile to remark that 
the order of magnitude of decoherence times in Eq. (|5l|) agrees with the one obtained 



by previous numerological arguments for gravity-induced localization [|l^]: " Although a 



detailed estimate of Tq would require a full theory of quantum gravity... it is reasonable 



to expect that for non-relativistic systems ..." What is new here in this regard is a 

fully defined dynamical model without any free parameter, which in principle allows for the 
explicit evaluation of any physically relevant quantity and for addressing crucial questions 
like the search for (gravitational-)decoherence free states of the physical operator algebra 



To be more specific, we have derived the first unified model for Newtonian gravity and 
gravity-induced decoherence. If the states \zj) in Eq. (^) are the pointer states of a 
measurement apparatus and \ej) are the measurement eigenstates of a microscopic system, 
the product state 

l^o) ® |ej) (53) 

i 

according to the traditional von Neumann model for the interaction between the two systems. 



is transformed into an entangled state 

J]c,|^,)®|e,). (54) 
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Obviously our previous analysis of the effect of the gravitational (self-)interaction on the 
quantum motion of the macroscopic body is not affected by the presence of the microscopic 
system, by which the reduction of the wave function occurs: 

CjCk \zj) \ej) (cfcl ® {zk\ — ^ ^ |cj|^ \zj) \ej) {ej\ ® {zj\ . (55) 

j,k j 



Of course one can look in principle for a collapse model [[18ip2| in terms of a stochastic 



dynamics for pure states, which, when averaged, leads to Eq. (|55D . Apart, in principle, 



from the non uniqueness of the stochastic realization |52|, stochastic models can certainly 



be useful as computational tools [53]. However the view advocated here considers density 
matrices arising from gravitational decoherence as the fundamental characterization of the 
system state and not just as a bookkeeping tool for statistical uncertainties. The fact that 
the apparent uniqueness of the measurement result seems to imply a real collapse is perhaps 
more an ontological than a physical problem, and presumably, if one likes it, that can be 
addressed by a variant of the many-world interpretation E 



VI. BLACK HOLE HEURISTIC 



Our first aim is to evaluate within our model the finite linear dimension of a collapsed 



matter lump, replacing the classical singularity. In order to do that we boldly use Eq. (|T5D 



for lengths smaller than /ig and /i2, namely in the limit /io,/i2 ~^ 0. This corresponds to 



the replacement of our meta-Hamiltonian with the model meta-Hamiltonian in Ref. p5 
where there is no gravitational interaction within observable and within hidden matter, 
while there is a Newton interaction between observable and hidden matter. This interaction 
is effective in lowering the gravitational energy of a matter lump as far as the localization 
length A ~ {H'^'E^ /GM^Y^^ is fairly smaller than the lump radius S. The highest possible 
density then corresponds roughly to (h^'E^ / (GM^)Y^^ = H , namely to 
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As to the space-time geometry, the Schwarzschild metric in ingoing Eddington- 
Finklestein coordinates {v,r,6,(j)) covers the two regions of the Kruskal maximal extension 
that are relevant to gravitational collapses p6|: 



ds^ = - [1 - 2MG/ (rc^)] dv^ + 2drdv + [dO^ + sin^ ed(j)'^] . (57) 
If in the region beyond the horizon we put x = v — Jdr[l- 2MG/ (rc2)]~\ then 

ds'^ =[l- 2MG/ (rc^)] dr^ - [l - 2MG/ (rc^)] dx^ + [dO'^ + sin^ ed(f)'^] (58) 



If we trust ( ^q) as the minimal length involved in the collapse, a future full theory 
of quantum gravity should include a mechanism avoiding the singularity at r = by the 
introduction of S as a regularization length. In particular, to characterize the region occupied 
by the collapsed lump, consider that for time-like geodesies at constant 6 and one can 
show that \dx/dr \ ~ r^/^ as r — > 0. This implies that the x coordinate difference Ax of two 
material points has a well defined limit as r 0, by which it is natural to assume that the x 
width of the collapsed matter lump is Ax ~ S. As to the apparent inconsistency of matter 
occupying just a finite Ax interval with d/dx being a Killing vector, one should expect on 
trans-Planckian scales substantial quantum corrections to the Einstein equations that the 
model gives on a classical level, with the dilaton and the ghost fields, though vanishing in 
the average, playing a crucial role. On the other hand we are proceeding according to the 
usual assumption, or fiction, of QM on the existence of a global time variable, at least in the 
region swept by the lump. In fact the most natural way to regularize ( pH]) is to consider it 
as an approximation for r > S of a regular metric, whose coefficients for r — correspond 



to the ones in (|58D with r = S, in which case there is no obstruction in extending the metric 
to r < 0, where taking constant coefficients makes d/dr a time-like Killing vector. As a 
consequence, the relevant space metric in the region swept by the collapsed lump is 

dslpACE ~ 2MG/ (Sc^) dx^ + -E? [de^ + sin^ ed(P'^] . 

The volume of the collapsed matter lump is then: 
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V ~ -E^Ax^MG/ (Hc2) = [n^l (GM^)] ^MG/c^ = ffM-'/ {G^c) . (59) 

According to the above view, thermodynamical equilibrium is reached, due to the gravita- 
tional interaction generating entanglement between the observable and hidden meta-matter, 
by which the matter state is a microcanonical ensemble corresponding to the energy 

E = Mc^ + GM^/E = Mc^ + GM^ [GM^/r?] ~ G'^M^/Jf, if M :^ Mp, (60) 

where Mp = ^JhcjG is the Planck mass, and to the energy density 

£ = EjV ~ G^cM^^/nJ. (61) 

For simplicity we treat the collapsed lump as a three-dimensional bulk, since treating it 
more properly, for the presence of the huge dilation factor in the x direction, as a string- 
like structure gives unchanged results. As this energy density corresponds to a very high 
temperature, not to be mistaken for the Hawking temperature, the matter can be represented 
by massless fields, whose equilibrium entropy is given by 

S ~ {Kb/ /v"^^'^) e^l^V = GM^Kb/ (he) . (62) 

Of course this result can be trusted at most for its order of magnitude, the uncertainty in the 
number of species being just one part of an unknown numerical factor. With this proviso, 
common to other approaches Eq. (^) agrees with B-H entropy. 

Our heuristic assumption of taking as gravitational energy of the collapsed lump just the 
expression given above in Eq. (|60|) is consistent with the connection existing between the 
temperature of the collapsed lump and Hawking temperature on purely thermodynamical 
grounds. In fact, if we take for granted that a future theory of quantum gravity will account 
for black hole evaporation, we can connect the temperature 

T ~ <fd^/KB ~ cGM^/KbH (63) 

of our collapsed matter lump with the (spectral) temperature of the radiation at infinity. If 
we model radiation by massless fields, emitted for simplicity at a constant temperature as 
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we are interested just in orders of magnitude, this temperature is defined in terms of the 
ratio Eoo/Soo of its energy Eoo and its entropy ^oo- It is natural to assume that, "once" 
thermodynamical equihbrium is reached due to the highly non-unitary dynamics close to the 
classical singularity, no entropy production occurs during evaporation, by which Soo = S. 
Then, if = Mc^ is the energy of the total Hawking radiation spread over a very large 
space volume, its temperature agrees with Hawking temperature, i.e. 

Too = (Eoo/E) T ~ {c'H/MGKb) . (64) 



VII. CONCLUDING REMARKS 



Of course the reversibility of the unitary meta-dynamics makes entropy decrease conceiv- 



able too , so that a derivation of the entropy-growth for a closed system (in principle the 
whole universe), in the present context, must have recourse to the choice of suitable initial 
conditions, like unentanglement between the observable and the hidden algebras. While the 
assumption of special initial conditions dates back to Boltzmann, only a non-unitary dynam- 
ics makes it a viable starting point, within a quantum context, for the microscopic derivation 
of the second law of thermodynamics, in terms of von Neumann entropy, for a genuinely 
closed system. This is meant without introducing generalized microcanocity conditions, and 



then renouncing isolation [5S 



It should be remarked that, for a realistic physical setting, most of the in principle 
observable degrees of freedom are yet out of our control and non-unitarity is the result 
of interactions with both fundamentally hidden degrees of freedom and with the environ- 
ment. Environment- induced decoherence [^], in most cases, may overshadow fundamental 



decoherence, even though the recent amazing experimental achievements in preserving and 
measuring quantum coherences make the detection of gravity-induced decoherence a less 
despairing task p6|j37| , |60| -p2| . In this respect the most natural experimental setting to look 



for gravitationally-induced decoherence seems to be that of Bose-Einstein (B-E) condensa- 



tion, due to the unprecedented scale of controlled quantum coherence achieved there ^3 
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In particular the localization mass-threshold is not too far from the present experimental 
limits, while its robustness with respect to the mass density variations may be a typical 
signature of our gravitational self-interaction. 

In conclusion, if we define a suitable non-unitary modification of fourth order gravity, by 
doubling the matter algebra and introducing a suitable constraint in order not to enlarge 
the observable algebra, we get the following outcomes: 

1) classical runaway solutions are absent and the ensuing classically stable theory may be 
made equivalent to Einstein gravity; 

2) the Newtonian limit is classically equivalent to ordinary Newton gravity; 

3) from a quantum viewpoint this non-unitary limit implies gravity induced localization 
and decoherence, which are compatible both with the wavelike behavior of microscopic 
particles and the classicality of the center of mass motion of macroscopic bodies; 

4) the model strongly supports the interpretation of the thermodynamic entropy of a 
closed system as von Neumann entropy and paves the way for the quantum foundations of 
the second law of thermodynamics; 

5) a bold use of the action at a distance limit of the model together with some geometric 
insight coming from Einstein gravity allows us to ascribe to the smoothed singularity of a 
black hole a finite entropy, which apart from an undetermined numerical factor coincides 
with the B-H entropy, and a very high temperature that is compatible with the much lower 
Hawking evaporation temperature. 
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